Domain enhanced interlayer coupling in ferroelectric/paraelectric superlattices 
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We investigate the ferroelectric phase transition and domain formation in a periodic superlattice 
consisting of alternate ferroelectric (FE) and paraelectric (PE) layers of nanometric thickness. We 
find that the polarization domains formed in the different FE layers can interact with each other 
via the PE layers. By coupling the electrostatic equations with those obtained by minimizing the 
Ginzburg-Landau functional we calculate the critical temperature of transition as a function of the 
FE/PE superlattice wavelength A and quantitatively explain the recent experimental observation of 
a thickness dependence of the ferroelectric transition temperature in KTaOa /KNbOa strained-layer 
superlattices. 

PACS numbers: 77.55.+f, 77.80.Dj, 77.80.Bh 



In the past decade refinements in deposition techniques 
have made it possible to fabricate nanoscale size oxide fer- 
roelectric superlattices with the objective to merge and 
optimize the technological properties of the constitutive 
materials Q, |^ . In designing such artificial structures 
an understanding of the physics of underlying processes is 
essential to determine whether the resulting characteris- 
tics are provided simply by the superposition of the bulk 
properties of the constituents or whether the interface 
and finite-size effects play a predominant role. 

Two competing types of phenomena that arise at the 
ferroelectric interface can affect the properties of the su- 
perlattices. The strain field, generated by the mechani- 
cal mismatch between the superlattice layers, infiuences 
the polarization orientation and generally increases the 
ferroelectric transition temperature Tc In contrast, 
the electric depolarization field, produced by interfacial 
surface charges is unfavorable to the formation of the 
ferroelectric phase In fact, in cubic perovskite-like 
ferroelectrics the situation can be even more complex 
due the formation of both 180° ferroelectric and 90° 
ferroelastic 0, Q, domains. Although the properties 
of ferroelectric superlattices can be governed by domain 
structure, no systematic study of this effect has to our 
knowledge been performed. 

In the present paper, we address the question of ferro- 
electric domain formation in a periodic superlattice struc- 
ture consisting of alternate ferroelectric (FE) and para- 
electric (PE) layers of equal nanometric width 2a f = 2ap. 
So as to avoid the complications of the effect of 90° fer- 
roelastic domains we assume that the ferroelectric layers 
have either natural or strain-induced c-oriented uniax- 
ial symmetry. We will show that the domain patterns 
formed in the different FE layers interact with each other 
across the PE layers via the spatially inhomogeneous de- 
polarization electric field emerging from the domains of 
the neighboring FE layers as shown in Fig. 1. This prox- 
imity type effect is dependent critically on the thickness 



i 


t 


I 


1 


i 


t 


i 


t 


w w v-* 




t 


f 


i 


! 

h4 


1 


f 


1 


t 



b) 



p 

T 



X 



FE 
PE 
FE 
>— PE 
FE 
PE 
FE 



Weak coupling 



Strong coupling 



FIG. 1: Domains and depolarization field in weak (a) and 
strong (b) coupled ferroelectric layers in FE/PE superlattice 
close to Tc 



of the PE layers. Our interest has also been motivated 
by a recent experimental study of FE/PE superlattices 
of KTaOs/KNbOa in which, as the superlattice wave- 
length A = 2a f + 2ap decreases, the ferroelectric tran- 
sition temperature Tc first decreases and then saturates 
below a certain layer thickness. Such dependence was re- 
produced using molecular dynamics simulations [loj but 
the role of depolarizing effects was not elucidated. At- 
tributing the observed behavior of critical temperature 
to the domain coupling in FE layers across the PE lay- 
ers, we calculate the dependence Tc{A) and show that 
this dependence correctly reproduces the experimental 
behavior. 

The physics underlying the domain formation and in- 
teraction is as follows. It is well known that a ferroelec- 
tric slab or thin film of thickness 2a/ will separate into 
domains. This is to reduce the energy of the depolar- 
ization field produced by the space and surface charges 
with charge density p{r) = divP that are provided by 
the discontinuity and non-uniformity of the polarization 
close to the crystal surface. In the Kittel approximation 
0,0,01, apphcable at T <C Tc, the polarization inside 
a domain is assumed to be uniform and almost equal to 
its equilibrium value Po- Thus the depolarization field is 
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proportional to the polarization discontinuity E^- ~ 47rPo 
and is confined to a thin layer of penetration length S that 
is roughly proportional to the domain size d. The equilib- 
rium domain structure is the result of a balance between 
two competing energies (calculated per unit surface of 
crystal): the depolarization energy F^ ^ SE^ ~ cIPq 
and the energy due to polarization gradients inside the 
domain walls F^ ~ a/APpn ~ APg(a//d) where the 
length scale parameter A is called a " domain wall width" 
and n ^ l/d is the wall concentration. Minimization of 
Fg + Fw gives the famous Kittel formula: d ^ -^Aa/ 

Qui 

The coupling between FE layers in a FE/PE superlat- 
tice is caused by the depolarization field emerging from 
the domain structure of the FE layers. This interaction 
is exponentially small and the FE layers are almost inde- 
pendent if the distance between them 2ap is larger then 
the penetration length 5. Taking into account that 5 



scales as 



we determine that for an equally layered 



superlattice with 2ap = 2a/ = A/2 this weak coupling 
regime is realized for long wavelength A. In the opposite 
limit of short A, the depolarization field penetrates the 
FE layers and so couples the domains in the neighboring 
ferroelectric layers. In this strongly coupled regime the 
domain size exceeds the superlattice wavelength and the 
superlattice behaves effectively as a uniform " composite" 
ferroelectric. 

Close to the ferroelectric transition the polarization 
varies gradually inside domains, the charge and depo- 
larization field penetrate FE layers so that the Kittel 
approximation is no longer applicable. In order to de- 
termine the domain structure parameters in FE/PE su- 
perlattices close to Tc and also the domain-formation- 
induced reduction of Tc with respect to the bare T^o of a 
bulk sample in which the depolarization field is screened 
by the short circuited electrodes, we generalize the anal- 
ogous calculations Ti'l for a thick ferroelectric plate. We 
start from the complete system of electrostatic equations: 



div(E + 47rP) = 0, rot E 0, 



(1) 



in which the polarization is related to the electric field by 
the constitutive relation P = P(E). This is to be deter- 
mined from both the nonlinear Ginzburg-Landau equa- 
tion for the z-component of the spontaneous polarization 

i f) 

in the uniaxial ferroelectric layers P — Pz : 
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(2) 



(here t = T/Tco - 1, ^ d^, + d^) and by the hnear 
relations for the two polarization components Px^} in the 
PE-layers as well as for the transversal component Px 
in the FE-layers: 



The susceptibility Sp of the PE-layers in (jj)) is assumed to 
be isotropic. The dimensionless parameter e|| ^ 1 in (0) 
is expressed via the Curie constant as e|| = C/Tco- The 
length ^0 — 6^ is estimated as the domain wall half- width 
at low temperatures, i.e. well below the phase transition. 
Hereafter we shall scale all the lengths in units of ^o- 

The formulation of the problem is completed by the 
electrostatic boundary conditions at the PE/FE inter- 
faces: 

i^(p)^„47r(Py)-Pi^)), Eif^^Ei^^ (4) 

and by the interface condition for the spontaneous polar- 
ization: 



dzP = AP, 



(5) 



where A is the extrapolation length P, that reflects 
the properties of the interface. 

Close to the ferroelectric transition Eq.© can be lin- 
earized by neglecting the term Pq^^P"^. The dimension- 
less transition temperature tc, where the instability to 
the formation of a multidomain structure first appears, 
can be found as the highest eigenvalue of the linearized 
system of equations , Q and with the appropriate 
boundary conditions l@J and (0). 

In terms of the electrostatic potential ip^-f'P^: Ei{i^^ = 
—dz.x'P^'^'^\ the full set of linearized equations takes the 
form: 

PE-layers: 



(9,2 + 9^)^(p)=0. 

FE-layers: 

^Trs^\t- V^)P = -dzip^f\ 
(^2 + ej.92)<^(/) ^ ^TrdzP. 

PE/FE interface: 

d,P = AP. 



(6) 

(7) 

(8) 

(9) 
(10) 



The eigenfunctions of the system of elliptic equations jnj, 
(|7|l and (|H1) are the linear superposition of the harmonic 
functions which we write as: 
For the upper (lower) FE-layer (see Fig. 1) 
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cos qx, 

(11) 
cos qx. 



For the central PE-layer 
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FIG. 2: Critical temperature tc — Tc/Tco — 1 (a) and domain 
structure wave vector q = n/da.tT = and T = Tc (b) as a 
function of superlattice wavelength A = 2af+2ap. Dotted line 
shows the asymptote I2H . Circles show the experimental data 
for KTaOa/KNbOs superlattice 0- The best fit parameters 
£|l = 400, £± = 500, Ep = 800, A = and A = 0.01 have been 
used. The length is scaled in units of ~ 6A 



In (|12() the boundary conditions H10|l were assumed. The 
periodicity of the oscillating factor cos qx reflects the for- 
mation of a regular domain structure in the x- direction 
with domain size d = ir/q. The solutions and H12|) 
can be periodically continued in the z-direction so as to 
follow the regular FE/PE superlattice structure. 

The parameters Pi, 2, fi,2 and ki_2 are found by sub- 
stitution of (inil and lO back into © , Q , I© , ® and 
(|10|l . This permits us to find two characteristic equa- 



J 



tions defining the eigentemperature t = t{q). 

Substitution of solutions lfTT)l into Eqs. Q) and © 
produces a homogeneous system of linear equations for 
(pi,2 and Pi, 2: 

47r£jj'"^ [t - kl ^ + 9^) Pi, 2 + fci,2 coth(fci,2a/) </?i,2 = 0, 
47rfci,2 Pi,2 - (fci,2 - £^q^) coth(fci,2a/) (^1,2 = 0, (13) 

that are compatible if the characteristic equation 

{kla - £±9") [t - kla + q^) + e\\kl^ = (14) 



is satisfied. 

Now, to obtain the t{q) dependence we need one more 
equation relating i, fc, and q. This equation follows from 
the boundary conditions. Condition 10 gives the follow- 
ing relation between (/3i,2 and Pi, 2: 

Qi cothfciO/ (^1 -I-Q2 coth/c2a/ V2 — 47r (Pi + P2) , (15) 
with Qi,2 — ki 2 + ep^cothgaptanh/ci ^ay, while condi- 
tion (O gives: 

Pi(fci tanhfcia/ - A) + P2(fc2 tanhfc2a/ - A) = 0. (16) 

Using H13|l we can express (pi,2 as a function of Pi, 2: 



</fl,2 



Pi 



47rfci,2 ^ , , 
1,2-^^2 „2 tanhfci,20/. 



'^1,2 - E-Lg 

Substitution of (|17|l into H15|l yields 
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kiQi 



kf - e±q'^ 



1 



k2Q2 



k2 — e±q'^ 



(17) 



0. (18) 



The compatibility criterion of equations H18|) and 1)16(1 
gives the second characteristic equation 



ki—e±q^ Ej^gtanhgOp -|- £p/c2 tanhfc2ay fc2tanhfc2aj 



fc| — ej^g^ s ±q tanh qup + Spkitanhkiaf fcitanhfeioy — A 



r 



(19) 



The eigentemperature t(q), transition temperature tc — 
maxg t(q) and the corresponding domain structure wave 
vector qc are found from equations 1(14(1 and ((19(1 after 
eliminating the variables fci2. The numerically obtained 
results for tc and qc as a function of superlattice wave- 
length A are plotted in Fig. 2 together with the exper- 
imental results of Ref. ^2] on tc{A) in KTaOa/KNbOa 
superlattice. The parameters we used to achieve this re- 



markably good fit were: 
A = and A = 0.01. 



e|l = 400, ej_ = 500, Sp = 800, 



weak and strong coupling limits are clearly seen in 
Fig. 2. With decreasing A the domain size decreases, 
goes through a minimum and then diverges. 

In the weak coupling regime (A > 20) the domain size 
is smaller then the superlattice wavelength and the de- 
polarization field is confined essentially at the FE/PE 
interfaces. Considering each EE layer as an independent 
thick film embedded in the PE media and using the sub- 
stitution p2 = ik2 (so that fc2tanhfc2a/ = — P2tanp2a/) 
we can simplify Eq. 1(191) to the form: 



Two regimes that correspond to the above described 



e^qtanhqap — SpP2 tanp2a/ ~ 



(20) 
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that was first worked out in Hi. Using Eqs. (O 
(|20|l the behavior of p2, (Zc ^^'^ reduces to: 



and 



P2 



tr. 




The critical temperature is inversely proportional to 
a/, exactly as was found numerically in 16]. 

In the strong coupling regime, below a certain critical 
thickness (A < 5) the domain structure abruptly dis- 
appears. The transition temperature of single-domain 
FE/PE superlattice is calculated from and ((T^ by 
assuming q — 0: 



a/' 



(22) 



or, equivalently, from the Landau energy of a periodic 
FE/PE structure with no depolarizing surface charges 
at the FE/PE interfaces. The FE/PE superlattice be- 
haves as a uniform composite ferroelectric with a critical 
temperature (|22|) , greater than that of the individual FE 
layers (|21|l . The positive surface A-term tends to increase 
tc and this is possibly the reason for the slight increase of 
the transition temperature in the KTaOs/KNbOa super- 
lattice at very small A 0, as shown in Fig. 2. Note that 
although the ferroelectric domains can exist in the strong 
coupling regime, their size is larger then the superlattice 
wavelength A and is defined by the global depolarization 
field of the sample. 

To estimate the temperature evolution of the domain 
structure we express its energy at T = (i.e. at t = — 1) 
in the Kittel approximation, assuming an abrupt struc- 
ture of the domain wall and a flat polarization profile 
inside the domains. The calculations are practically the 
same as those performed for the domain structure of the 
uniaxial ferroelectric film surrounded by paraelectric pas- 
sive layers and embedded in a short circuited capacitor 
^ . The free energy is the sum of the domain walls energy 
and the electrostatic contributions: 



ttF/P^ = 2afAq - 
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ra=l,3,. 
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where 



9niq) 



1/2 



(23) 



9niq) — EpCothnqap + ^eiiex coth — I nqaf 

\^\\ 

and the "domain wall width" A can be found by in- 
tegration of the Landau energy of the wall as A = 
ifi /^(tanh''(x/V2) - l)dx « 87rV2/3e||. This takes 

into account the fact that P{x) = Pg tanh(a;/\/2) is the 
exact single- wall solution of Eq. ||2J) at < = — 1. The re- 
sult of the numerical minimization of 1231) is given by the 
dashed hue in Fig. 2b. When gcOp, {e^/e^^Y^^qcaf » 1 
it is approximated by the generalized Kittel formula: 



21C(3) 1 



21C(3) 

Ep + y/EfET 27rV2 2a/ 27r^/2 



Since, as shown in Fig. 2, the plots q{A) at T = and 
aXT = Tc practically coincide, we conclude that the tem- 
perature dependence of the domain structure wave vec- 
tor is very weak. The complete calculation of the domain 
structure evolution and of its dielectric constant over the 
entire temperature region will be published elsewhere. 

To conclude, we have demonstrated that uniaxial fer- 
roelectric domains can substantially influence the prop- 
erties of the FE/PE superlattices. Depending on the 
wavelength A, the superlattice can be in different domain 
states. For large A (typically > 5 — 15nm) each FE layer 
has an independent domain structure. At smaller A the 
domains in neighboring FE layers interact through the 
PE layers via the emerging depolarization field and this 
results in a dramatic increase of the domain width. In 
this regime the superlattice structure behaves as an effec- 
tive composite uniform ferroelectric where the large-scale 
domains penetrate throughout the entire sample and are 
governed by the global depolarization field. 

We have calculated the ferroelectric transition temper- 
ature as function of A and have explained the recently ob- 
served small- A saturation of Tc in KTaOs/KNbOs super- 
lattice by crossover to the regime of strongly coupled FE 
layers. We have also shown that the critical superlattice 
crossover wavelength is nearly temperature independent. 
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